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Abstract 

Background: The spectroscopic conductivity distribution of tissue can help to explain 
physiological and pathological status. Dual frequency conductivity innaging by 
connbining Magnetic Resonance Electrical Property Tomography (MREPT) and 
Magnetic Resonance Electrical Impedance Tomography (MREIT) has been recently 
proposed. MREIT can provide internal conductivity distributions at low frequency 
(below 1 kHz) induced by an external injecting current. While MREPT can provide 
conductivity at the Larmor frequency related to the strength of the magnetic field. 
Despite this potential to describe the membrane properties using spectral information, 
MREPT and MREIT techniques currently suffer from weak signals and noise 
amplification as they both reply on differentiation of measured phase data. 

Methods: We proposed a method to optimize the measured phase signal by finding 
weighting factors according to the echo signal for MREPT and MREIT using the ICNE 
(Injected current nonlinear encoding) multi-echo pulse sequence. Our target weights 
are chosen to minimize the measured noise. The noise standard deviations were 
precisely analyzed for the optimally weighted magnetic flux density and the phase 
term of the positive-rotating magnetic field. To enhance the quality of dual-frequency 
conductivity images, we applied the denoising method based on the 
reaction-diffusion equation with the estimated noise standard deviations. A real 
experiment was performed with a hollow cylindrical object made of thin insulating film 
with holes to control the apparent conductivity using ion mobility and an agarose gel 
cylinder wrapped in an insulating film without holes to show different spectroscopic 
conductivities. 

Results: The ability to image different conductivity characteristics in MREPT and MREIT 
from a single MR scan was shown by including the two objects with different 
spectroscopic conductivities. Using the six echo signals, we computed the optimized 
weighting factors for each echo. The qualities of conductivity images for MREPT and 
MREIT were improved by optimization of the phase map. The proposed method 
effectively reduced the random noise artifacts for both MREIT and MREPT. 

Conclusion: We enhanced the dual conductivity images using the optimally 
weighted magnetic flux density and the phase term of positive-rotating magnetic field 
based on the analysis of the noise standard deviations and applying the optimization 
and denoising methods. 
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Background 

The conductivity spectra of biological tissues can provide diagnostic medical information 
from the estimation of physiological and pathological conditions of in-vivo and ex-vivo 
tissue. However it is difficult to produce high resolution conductivity images inside the 
human body [1,2]. The conventional conductivity imaging methods have limited spatial 
resolution and sensitivity inherited from the ill-posed nature of the problem [3]. In order 
to achieve sensitive conductivity images with high resolution, electric impedance imag- 
ing techniques based on a magnetic resonance imaging (MRI) have been propsed. These 
include magnetic resonance electrical impedance tomography (MREIT) and magnetic 
resonance electrical property tomography (MREPT) which are under active investigation 
[4-17]. Both methods use the internal magnetic field information obtained from the phase 
data of MRI scanner to reconstruct the internal conductivity image. 

MREPT can provide the electrical conductivity information at the Larmor frequency 
by measuring the phase of positive rotating field due to the applied Bi field. It does not 
require application of external current and therefore directly recovers the conductivity 
distribution by taking the derivative of the measured phase signals twice. MREIT needs a 
pair of electrodes to inject current into an imaging object during MRI scan for measuring 
a magnetic flux density induced by the external injecting current. The MREIT technique 
used only the ;2-component of magnetic flux density, Bz, of B = {Bxy By, Bz) to reconstruct 
the cross-sectional apparent conductivity image at a lower frequency range (below 1 kHz) 
[18-22]. The phase difference approach with an interleaved encoding scheme was adopted 
to cancel the systematic artifacts accumulated in phase signals and also reduce the ran- 
dom noise artifacts. Recently, a simultaneous conductivity imaging technique using a 
combination of MREPT and MREIT was proposed to provide the dual-frequency con- 
ductivities of tissue from a common MR scan [23]. MREPT using -mapping technique 
visualizes the conductivity and permittivity distributions at the Larmor frequency and 
MREIT recovers the apparent conductivity distribution when injecting low-frequency 
external current through the attached electrodes. Since the biological tissues show the 
frequency dependent conductivity property [1,2], the simultaneous dual-frequency con- 
ductivity imaging using a single MR scan is beneficial to provide distinct electrical features 
of tissues quantitatively. 

Both conductivity imaging techniques use the phase signals of measured MR data. 
MREIT and MREPT commonly suffer from weak signals and noise amplification by using 
the derivative of the measured signal. The noise level of B^ in MREIT is inversely pro- 
portional to the signal-to-noise ratio (T^') of the MR magnitude image and the current 
injection pulse width. Because of the small amount of injection current and the poor qual- 
ity of measured Bz change due to the injected current, it is difficult to perform in vivo 
human experiments using a conventional MR pulse sequence. To enhance the magnetic 
flux density due to the injected current in MREIT, the injected current nonlinear encod- 
ing (ICNE) method was introduced. This extended the duration of injecting current until 
the end of a readout gradient, and improved the signal by using a multi-echo train MR 
pulse sequence [24-26]. MREPT also suffers from low sensitivity due to the inherently 
poor signal to noise ratio and noise sensitive characteristics as it also needs the derivative 
of measured data, is very sensitive to the measured noise. 

In this paper, we adopt a multiple spin echo MREIT pulse sequence based on the ICNE 
scheme to measure multiple phase data for MREIT and MREPT images. The acquired 
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multiple phase data can be decomposed to the phase term reflecting the magnetic flux 
density signal induced by the injected current and the other phase term of positive rotat- 
ing field due to the applied Bi field. We analyze the noise level of the two decomposed 
phase terms and minimize the measured random noise artifacts by applying the optimal 
combination of multiple phase terms. Also, we apply a denoising technique to the opti- 
mized magnetic flux density for MREIT and to the phase signal for MREPT in order to 
improve the quality of the reconstructed conductivity images. We prepared a conduc- 
tivity phantom consisted of two different kinds of anomalies to show the difference of 
MREIT and MREPT. A phantom experiment is conducted to validate that the proposed 
method is able to improve the qualities of reconstructed dual-frequency conductivity 
images compared to the results of using conventional MREIT and MREPT reconstruction 
algorithms. 

Methods 

Governing equation 

We denote an imaging object as ^2. The admittivity in ^2 is /c = a {oj) + iooe (oj), where a (o)) 
and 6 {o)) are the conductivity and permittivity, respectively, at the angular frequency co. 
Time-harmonic Maxwell's equations relate the electric field E and magnetic flux field H: 

V X E(r) = -/a;/xoH(r) and V x H(r) = /c(r)E(r), r G ^2, (1) 

where the current J and the electric field intensity E satisfy the relation J = /cE by Ohms 
law. 

In MREIT, we can assume that the injecting current is sufficiently low frequency to 
meet i(0€{(o) ^ 0 and neglect icoiJLo ^ 0, Maxwell's equation (1) satisfies 

-V^HW = V X V X H(r) = V x (ai(r)E(r)) 

= Vori(r) X E(r) = x (V x H(r)), r e Q. 

where ai = a(co) denotes the low-frequency electrical conductivity responding to the 
externally injecting current. Since the electric field intensity E = — Vw is a gradient form 
under negligible angular frequency co, we have 

^ V • (ai^uir)) = 0 in Q 

(3) 

—aL^u(r) ' n(r) = ^(r) ondQ 

where n is the outward normal vector on the surface dQ and ^ is the current density on 
the surface. 

In MREPT, by assuming the locally homogeneity, the magnetic field H can be expressed 
as a simple form: 

-V^HCr) = V/c(r) x E(r) - /66>/xo/c(r)H(r) = -/66>/xo/c(r)H(r), r e Q. (4) 
The phase terms retrieved from the measured phase are the combination of the phase 

^ -(-/// j-[ iJ-[ 

terms of positive //+ = ^ ^ negative H~ = ^ ^ rotating field. For a restricted 
situation such as the usage of transmit-receive coil, the phase term of positive rotating 
field can be retrieved from the measured phase. The transverse field of H can be decom- 

H -\-iH 

posed into the positively rotating field /f+ = 2 ^ negatively rotating field 
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H- = We denote (p+ and (p as the phase term of and H , respectively. When 

single channel transmit-receive coils are used, the transmit can be reasonably similar 
to (p~ [16]. The measurable positively rotating field = ^""^^^y of H = (Hx,Hy,Hz) 
also satisfies the relation under the locally homogeneity condition 

-V^//+(r) = -icoiioK{r)H^{r), r e (5) 

By separating the real and imaginary parts of k and assuming (p~^ ^ (p~, the phase term of 
//+ can estimate the conductivity gh = cf {o)) which reflects the high-frequency property 
responding to the RF pulse [16]: 

OHix) ^ ^^lAZ, r G a. (6) 
ICNE multi-echo pulse sequence and multiple phase data 

The ICNE multi-echo pulse sequence is adopted to acquire multiple echoes measured at 
the echo time Tej, j = 1, • • • , A/^, where Ne is the number of echoes per repetition time 
Tr, Figure 1 shows its schematic diagram where and /~ are sequentially injected pos- 
itive and negative currents synchronized with the applied 180° rephasing pulses with Ne 
times. The ;-th complex image, corresponding to the ;-th readout gradient generated 
with the positive and negative injecting currents and /~, can be expressed as 

i;f{v) = M«^,/r)e'^^e'^"e±(-i>''>^c/5./W, / = i, . . . ,Ne (7) 

where M^^y(r) = M^y(r)e ^^/^^^''^ is the y-th transverse magnetization at a flip angle 
Of, 8^ is the systematic phase artifact, y is the gyromagnetic ratio of hydrogen, T2 is the 
transverse relaxation time, and Tj is the ;-th duration of the injecting current. 
The y-th transverse magnetization M^^y(r) can be rewritten in detail as 

M^^/r) = CiMo(r)//-(r) sin(C2a|//+(r)|)e"^, ; = 1, • • • , Afe (8) 

where Q and C2 are system-dependent constants. Assuming 8^ ^ 0, the phase signal can 
be separated as 

Vf(r) := arg (^^±(r)) = 2^+(r) ± {-ly yT,jB,j{v\ j=h'".NE. (9) 
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Figure 1 Schematic diagram for ICNE multi-echo pulse sequence to acquire dual-frequency 
conductivity images. 'RF' shows the timing for excitation of 90° and 1 80° RF pulses. Gz, Gy and Gx are slice 
selection gradient, phase encoding gradient and readout gradient, respectively. is the y-th echo signal. 
'ICNE' shows the alternating injection currents in the form of pulses synchronized with the RF pulse. 
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By subtracting and adding the phase signals Vp, we can obtain the magnetic flux density 
Bz induced by the externally injecting current at the low-frequency and the phase term of 
at the Larmor frequency simultaneously. 

VUv) + 7^r(r) = 4^+(r), ; = 1, • • • .Ne 

\ ' ' (10) 

Vf (r) - Vr (r) = 2y T.iB^i (r), ; = 1, • • • , iV^ 

From the relation 10, we can compute the /-th phase signals for MREPT and the /-th 
;^-component of magnetic flux density for MREIT: 

<pI = ^^, j=l,---.NE 

(11) 

V--V~ 

The noise standard deviation of the measured magnetic flux density in 11 is inversely 
proportional to the current injection time and the signal-to-noise ratio (T^') of MR 
magnitude image following [27,28] 

sdBj(r) = ^— , ; = 1,...,7V£. (12) 

2yr/T/(r) 

A similar analysis provides the noise standard deviation of (pj^ as 

sd^+(r) = — -. — , ; = l,...,Afe. (13) 
Optimal combination of multiple B^j 

In order to enhance the conductivity image at the low frequency using MREIT, we have to 
improve the quality of the Bz signal An alternative to reduce the noise of Bz is to find an 
optimal combination of multiple B^j, y = 1, • • • , Ne> by determining the weighting factors 
which satisfied conditions of fy(r) > 0 and Yl^i ?/ (^) = 1- When we assume that the 
measured /c-space signals are contaminated with independent identically distributed (IID) 
complex Gaussian random noise as same amount of expected noise standard deviation of 
each B^j in 12, the noise variance of the weighted B^j satisfies the following relation 

Var(4(r)) = ^ ^^(r)Var(5,;(r)) a ^ , /.^ . (14) 
j=l j=l ^J^V 

To determine the optimized weighting factor §y(r), the noise variance of the combined B^ 
is required to minimize the following conditions: 

^H=1.-,Ne ^^^^ 
subject to GB.i^i. • • • , ^Ne)(^) = h §;(r) > 0 

where := Effi and > > > := E^i^/W. The 

method of Lagrange multipliers to solve the optimization problem in 15 introduces a new 
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variable X and requires the condition that the gradients of Gbz parallel Using 

the method of Lagrange multipliers, the optimal weighting factor §y(r) can be obtained as 

where vI/^-(r) := T^. |^±(r)p. 

The optimized magnetic flux density to produce a conductivity image in MREIT can 
be generated by the weighted average of multiple B^j with the determined weighting fac- 
tor of §y(r) as Bl(r) = Y^=i ^j(^)^z^(^) using the ICNE multi-echo pulse sequence. The 
noise standard deviation of can be precisely described by substituting 16 into the noise 
variance of bI in 14: 

— ^ss* ■ 

where AfMxy denotes the noise level of magnitude image. 
Optimal combination of multiple (p^ 

The noise variance of the weighted cp^ satisfies the following relation similar as 14. 

Ne Ne ^2^J,^^ 

Var(^+'^) = J2 x/(r)Var(^+) oc ^ -^—^ (18) 
j=l j=l 

The noise variance of the combined cp^ with weighting factors is required to minimize 
the following conditions, 

min J'^+(xi.--- .XAfe)(r) 
x./=i,-,A^. (19) 

subject to Q^p+ (xi, • • • , XNe)(x) = h X;(r) > 0 

where J'^+(xi, • • • , XAfe) := E/ii j^^^ and G^+ixiy • • • , XAfe) := E/ii X;(r). The opti- 
mal weighting factor, x/(i')» corresponding to each echo phase can be determined by using 
the method of the Lagrange multipliers as 

««=-|f- (20) 

where ^y(r) := \^j^(r)\'^. 

The optimally weighted average of multiple phase signal cp^ with the determined 
weighting factor of Xj(^)> ^"^'^(r) = E/^i Xj(x)^f(x)y can produce an enhanced con- 
ductivity image in MREPT. The noise standard deviation of the optimized phase signal, 
cp^'^y is computed by substituting (20) into the noise variance of ^"^'^ in (18) to estimate 
pixel-by-pixel noise precisely: 

sV^(r) = ^^2L_ (21) 

Considering only random noise effects, the optimal weighting factor Xj(x) reduces the 
noise level, depending on the echo number and 72 -decay rate at each imaging pixel. 
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Denoising method using the estimated scf^i and scf^+,x 

A denoising technique is applied to gradients of the optimized magnetic flux density for 
MREIT image and the optimized phase signal for MREPT image since the reconstruc- 
tion procedures only require the differentiated measured data. Also, it is advantageous to 
denoise V^f and Wcp^'^ because the measured magnetic flux density of and the phase 
of are inherently continuous without conventional edge information. 

To denoise the optimized V^f and V^+'^, we use the following conventional reaction- 
diffusion equation [29,30]: 

^ ^ (22) 

v(r,0) =/(r) 

where the parameter ^ is a fidelity term and the function \l/ satisfies 

if(s)>0, hmxlf(s) = 0. (23) 

Here, the initial state / denotes each component of the optimized magnetic flux density 
VbI or the phase signal V^+'^, which is to be denoised. The time dependent solution 
v(r, t) is the denoised image of the initial state image / 

To remove the random noise artifact while preserving the edge information, it is impor- 
tant to determine the diffusion function ^ and fi. Using the estimated noise level of 5/i^+,x , 
the diffusion function i// and the fidelity function can be specified as 

xlf(r) (X \ and fi(r) oc sd^+,x (r) (24) 

|Vv(r,/:)l + ^6 

where is a small parameter to guarantee the positive sign of |Vv(r, 01 + ^e- For the 
magnetic flux density ^f, the diffusion function ^1/ can be selected similarly. 

The quality of magnetic flux density signal in MREIT depends on the magnitude inten- 
sity and width of injected current simultaneously, while the quality of phase signal in 
MREPT only depends on the magnitude intensity. Therefore, there are common effects 
due to the decay of the magnitude intensity in MREPT and MREIT processing. There are 
also different effects caused by externally injecting current in MREIT. These occur when 
we calculate each of the optimal weighting factors and the noise standard deviation in 
the optimal combination of multiple phase and magnetic flux density images. The pro- 
posed optimization method uses the decay rate of magnitude intensities to determine the 
weighting factors in MREIT and MREPT, which relate with the T2 values. However, the 
determined weighting factors in 16 and 20 only include the measured magnitude inten- 
sity at each echo time T^^., / = 1, • • • , AT^. The proposed method therefore does not need 
estimates T2 values to optimize the multiple phase signals. 

Phantom design and experimental setup 

In order to evaluate the proposed optimizing method using ICNE multi-echo MR pulse 
sequence, we designed a phantom which can show different characteristics of recon- 
structed conductivity images by MREIT and MREPT. Figure 2(a) and (b) show the 
phantom configuration at the middle slice of the object and the photo view of the 
phantom, respectively. We built a cylindrical phantom filled with a saline of 0.2 Sm~^ con- 
ductivity (0.3 g/L NaCl and 1 g/L CUSO4). The diameter and height of the phantom were 
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Figure 2 Experimental set-up for a phantom, (a) The phantom configuration at the middle slice of the 
testing object, (b) a cylindrical phantom with saline solution of 0.2 Sm~^ including a thin film object and an 
agarose anomaly wrapped by an insulating film, and (c) a magnitude image at the middle slice acquired at 
TE = ^S ms. 



11 cm and 14 cm, respectively. We attached four carbon hydrogel electrodes (HUREV Co. 
Ltd, Korea) on the side of the acrylic container for injecting current. 

Two different objects were positioned inside phantom. The left one was a thin hollow 
cylindrical object with the diameter of 4 cm using an insulating thin film of 0.4 mm thick- 
ness. We punched four holes which had 2 mm diameter with equally spaced around the 
circumference. Because we filled the same saline of 0.2 Sm~^ inside and outside of the 
hollow cylindrical object with holes, the apparent electrical conductivities inside and out- 
side of object were determined by the movement of ions through the holes [31]. The right 
cylindrical object with the diameter of 3 cm was made of an agarose gel (1 g/L CUSO4, 
2.1 g/L NaCl, 15 g/L Agar) to generate different conductivity, spin density and 72 -decay 
compared with the background saline. The conductivity of agarose gel was 1.10 Sm~^, It 
was wrapped in a thin insulating film without holes. 

A transversal injecting current of 10 mA was introduced into the phantom via a pair of 
recessed carbon hydrogel electrodes attached at the middle of the phantom. Figure 2(c) 
shows the magnitude image at the middle slice. We used the ICNE multi-echo MR pulse 
sequence by injecting current of alternating polarity synchronized with the multiple 
refocusing pulses as shown in Figure 1. 

Imaging parameters in a 3T MRI scanner (Achieva, Philips) with birdcage transmit- 
receive (Tx/Rx) RF head coil were as follows: repetition time Tr = 1200 ms, data 
acquisition time width Ts = 3.584 ms, echo-spacing Te^^ = 15 ms, total number of echo 
Me = 6 and number of averaging = 4. The reconstructed image matrix was 128 x 128, 
with a FOV of 180 x 180 mm^, 5 mm slice thickness. Total scanning time for both verti- 
cal and horizontal injection currents was 20 minutes with an interleaved phase encoding 
acquisition. 

Results 

Figure 3 shows the multiple MR magnitude images measured at the echo time Tej = 
15 X y ms for ; = 1, • • • ,6, where the right anomaly shows rapid decay of the transverse 
magnetization due to magnetic field non-uniformity and spin-spin transverse relaxation. 

Figure 4 shows the phase images of of //+ map and magnetic flux density B^j 
induced by the vertically injected current using ICNE multi-echo pulse sequence at the 
echo time Tej = 15 xy ms,y = 1, • • • , 6. The images in Figure 4(a) and (b) were estimated 
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40 



Figure 3 MR magnitude images measured at each echo time. Echo number marked on the upper left 
part of each image as (1), (2), • • • , (6). Magnitude images using ICNE multi-echo pulse sequence acquired at 
the echo time Tej = 15 x ; ms for; = 1 , • • • , 6. 



by adding and subtracting the phase signals of ^j^, respectively. Since the electromagnetic 
wave at the Larmor frequency of 128 MHz could penetrate the thin insulating film and 
there was same saline inside and outside of the thin insulating wall of the left hollow cylin- 
der, the measured phase values of cp^ for the left object in MREPT could not produce 
conductivity difference between inside and outside. However, the measured B^j images at 
the low frequency in MREIT showed the different value inside the left anomaly because 
MREIT measures the apparent magnetic flux density corresponding to the ion mobility 
and intrinsic conductivities of composite materials. Note that the slope of magnetic flux 
density was changed abruptly due to the high current density around the holes in the same 
direction of the current injection. Regarding the different conductivity agarose gel cylin- 
der wrapped in the thin insulating film without holes, both images in Figure 4(a) and (b) 
could present distinctions between inside of the object and the background saline. 

Using the proposed optimization method with the estimated noise standard deviations 
sd + and sd j = 1, • • • ,6, we found the optimal weighting factors for the phase of 
and the magnetic flux density by the injected current, respectively. Figure 5(a) and 
(b) show the determined weighting factors for dual-frequency conductivity images in 
MREPT and MREIT combined method, respectively. In MREPT, the intensity of weight- 
ing factors was monotonically decreased because the noise level of cp^ only depended on 
r2 -decay rate of ^j^, j = 1, • • • ,6, whereas the weighting factors in MREIT show the dif- 
ferent characteristics because they depended on the combination of 72 -decay rate and 
the duration of injected current simultaneously. In Figure 5(b), the intensity of weighting 
factors in the background area was monotonically increased because the effect of time 
increasing for the injection current was greater than the magnitude attenuation effect of 




(b) 




Figure 4 Phase images and magnetic flux density images of each echo, (a) Phase images cpj^ of /-/+ map 
in radian and (b) magnetic flux density images B^j in 'nT' unit using ICNE multi-echo pulse sequence at the 
echo time If. = 1 5 x ; ms for; = 1 , • • • , 6. 
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0 

nO.4 



1 

Figure 5 Images of weighting factors corresponding to each echo, (a) Images of weighting factors 
corresponding to eacli eclioes for (p'^' and (b) images for magnetic flux density B^j to produce an optimal 

<^+'^ and b\, respectively. The weighting factors were scaled from 0 to 1 . 



echo signals. However, the weighting factors in the right agarose gel anomaly region were 
different due to the combination of both effects. The weighting factor values at the second 
echo were higher than those of other echoes. 

Figure 6 shows the reconstructed dual-frequency conductivity images using the multi- 
ple echoes (p'^ and B^^^ for j = 1, • • • ,6 and d = 1,2, where B^j and B^j were the measured 
magnetic flux density data by the vertically and horizontally injected currents, respec- 
tively. To reconstruct the high-frequency conductivity distribution in Figure 6(a), we used 
phase for each echo signal as 

V2^+(r) 

cr^H(r) = ^ , rea ; = 1,...,6. (25) 

Here, we could not detect the hollow cylinder in the left side made of the thin insulating 
film with holes. The high-frequency conductivity inside and outside of the left object was 
estimated to be the same value. The quality of reconstructed conductivity image for each 
echo in Figure 6(a) was monotonically degraded corresponding to the weighting factors 
shown in Figure 5(a). Table 1 shows the values of weighting factor inside and outside of 
the right agarose object wrapped by insulating film. 





S/m . 


(a) 


^^^^^^^^^^ ""^^1 


1 -n ^ 




S/m 


(b) 


(1) (2) (3) (4) (5) (6) 

5^ -f -f 




Figure 6 Reconstructed conductivity images of each echo, (a) Reconstructed conductivity images using 
the multiple echoes for; = 1 , • • • , 6 at Larmor frequency and (b) conductivity images using 8^^- for 
j = 1 , . . . , 6 at the low frequency. 
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Table 1 Weighting factor values for MREPT inside (7^/») and outside (IZout) the agarose 
anomaly 

Ne = 1 Ne = 2 Ne = 3 Ne = ^ Ne = 5 Ne = 6 

IZout 0.2030 0.2046 0.1737 0.1570 0.1385 0.1232 

Hin 0.4926 0.2830 0.1215 0.0617 0.0268 0.0143 



We estimated the projected current density J^, which was the optimal divergence-free 
component from the measured ;^-component magnetic flux density in order to recon- 
struct the low-frequency conductivity images [32]. The projected current density consists 
of the background current density = — Vof and V^yy = {^y—^^^y where the 
potential r) satisfies at the imaging slice ^2^ C ^2 



xy ^ 



= 0 



on dQf 



(26) 



where Wxy denotes the two-dimensional gradient. The absolute conductivity distribution 
can be written as a matrix form using the two measured Bf, d = 1,2, data 



(Ao + A^ 



^ \dx ay / 



in Qt 



(27) 



where is an imaging slice, r = log a, and 



Ao = 



dy 
dy 



dx 

dx / 



A* = 



dx 
\ dx 



dy 

drf^ 
- dy ) 




(28) 



Figure 6(b) shows the reconstructed low-frequency conductivity images by solving 27 
using the /-th measured two magnetic flux densities 5^^, d = \,2 obtained with the ver- 
tically and horizontally injected currents. The reconstructed conductivity values inside 
the hollow cylinder with holes in the left side reflected the externally injected current 
flowed through the four holes. Regarding the agarose gel object wrapped in the insulat- 
ing film without holes, the apparent conductivity reconstructed by solving 27 shows the 
same anomaly as an insulator since the measured Ez data only reflected the current den- 
sity distribution due to the externally injected current through the attached electrodes. 
The quality of conductivity image using the sixth echo data was poor inside the agarose 
gel anomaly region due to relatively short 72 -decay relaxation time. It was consistent with 
the corresponding weighting factor in Figure 5(b). Table 2 shows the values of weighting 
factor inside and outside of the right agarose object wrapped by insulating film. 

Figure 7(a) and (e) show the normalized noise standard deviations of optimally weighted 
ip'^'^ and ^1, respectively. In previous approaches without using multi-echo data with 
optimizing method, the conductivity images for MREPT and MREIT were reconstructed 
as shown in Figure 7(b) and (f). Due to the different 72 -decay rate and the duration of 



Table 2 Weighting factor values for MREIT inside (7^/») and outside {IZout) the agarose 
anomaly 





Ne = 1 


Ne = 2 


Ne = 3 


iV£ = 4 


Ne = 5 


Ne = 6 


^out 


0.0238 


0.0734 


0.1269 


0.1936 
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0.3239 




0.1404 
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0.2158 


0.1849 


0.1214 


0.0910 
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Figure 7 Optimized dual-frequency conductivity images with and without applying the optimization 
and denoising method, (a) and (e) Normalized noise standard deviation of tlie optimally weighted phase 
of /-/+ and fif, respectively, (b) and (f) Reconstructed conductivity images without using the optimization 
and denoising method, (c) and (d) Reconstructed conductivity distributions at Larmor frequency using the 
optimized (p^'^ without denoising and with the proposed denoising technique applying to V<^+'^, 
respectively, (g) and (h) Reconstructed low-frequency conductivity distributions without denoising and with 
the proposed denoising technique applying to VB^'^ , d = 1,2, respectively. 



injected current, the intensities of sd^p+a and sd^^ represented different characteristics. 
The reconstructed dual-frequency conductivity images with the optimally weighted q)^'^ 
and , d = 1,2, were displayed in Figure 7(c) and (g). Solving the reaction-diffusion 
equation 22 with the estimated noise standard deviations sd(p+,x and sd^^ , we denoised 

V(^+'^ and V^f'^; the iteration number was 400 and the regularization function of P was 
j\f(^d\ X ) MREPT case, where J\f{sd^p+,x) was the normalized noise standard devi- 
ation in Figure 7(a). The denoising parameters for MREIT were the same as those for 
MREPT. Figure 7(d) and (h) show the reconstructed dual-frequency conductivity images 
after denoising of V(^+'^ and V^f'^, respectively. 

Discussion 

MREPT and MREIT used an internal magnetic flux density information from a phase 
imaging can provide conductivity images with high resolution compared to the conven- 
tional electrical impedance imaging techniques. Moreover, the interleaved measurement 
technique with alternating injected current polarity based on a multi-spin-echo pulse 
sequence, enables simultaneous dual conductivity images. From a single MR scan, 
we could obtain the low-frequency conductivity images by subtracting and the high- 
frequency conductivity images by adding the measured phase terms. Despite their 
advantages and usefulness, MREPT and MREIT techniques commonly suffer from noise 
amplifications in phase imaging. The proposed method used a ICNE multi-echo pulse 
sequence based on a spin echo pulse sequence in order to suppress the background field 
inhomogeneity. This also influenced to the reconstructed high-frequency conductivity 
using the MREPT technique. Even though the estimated weighting factors for MREIT and 
MREPT optimally reduce the random noise artifacts, non-uniform noise artifacts from 
various sources can severely deteriorate the reconstructed conductivity distributions. To 
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be a clinical device used practically, it is important to develop noise reduction techniques 
taking into account the physical properties and experimental environments of MREIT 
and MREPT. 

The designed phantom consists of two different kinds of anomalies to show the fre- 
quency dependent conductivity information. The stable apparent conductivity contrast 
of hollow cylindrical object made of thin insulating transparency film was only controlled 
by ion mobility through holes on the film wall. It excluded other effects caused by any ion 
concentration gradient. Therefore, MREIT images could only present the hollow cylin- 
drical object with holes when comparing the phase terms and V~ in 10. The apparent 
conductivity of the region inside cylinder was a non-zero conductivity value because the 
externally injected current entered the anomaly through the holes. In MREPT, the recon- 
structed high-frequency conductivity was same inside and outside of the anomaly since 
the measured phase of only reflected the material conductivities inside and outside 
of the cylindrical film. Another agarose gel cylinder wrapped in an insulating film with- 
out holes shows different conductivity characteristics in MREPT and MREIT. Due to the 
short 72 -decay relaxation time of agarose gel, the noise levels of 5(i^+,x and sd^^ were 
relatively severe compared to those in the other regions. The noise level of each sd + 
was strictly increasing with respect to the echo number depending on the 72 -decay rate, 
whereas those of sds^ had different characteristics depending on T^j and the 72 -decay 
rate simultaneously. 

MREIT provides the conductivity distribution at a low frequency, whereas MREPT pro- 
duces the conductivity at the Larmor frequency of 128 MHz at 3 T. The cell membranes 
consisted of the phospholipid bilayer with embedded proteins behave as a capacitor or 
insulating film. It shows a complicated pattern of conductivity depending on the mea- 
sured frequencies. For the biological tissues, the membranes restrict the flow of current of 
low frequency. MREIT images reflect the membrane properties quantitatively and MREIT 
has a potential to visualize the anisotropic conductivity tensor map. On the other hand, 
the membranes become transparent at high frequency and provide a relatively degraded 
sensitivity in the conductivity image. MREPT may provide different conductivity charac- 
teristics compared to MREIT results. Based on the previous studies related to the complex 
conductivity spectra measured at multiple frequencies within the range of 10 Hz to sev- 
eral kHz, spectroscopic complex conductivity distribution can contribute to explain the 
physiological and pathological status of internal tissues [33]. The dual-frequency conduc- 
tivity imaging by combination of MREPT and MREIT is an advanced technique to show 
different information in spite of displaying two extreme cases at low and the Larmor fre- 
quencies. It is meaningful to develop a method to produce multi-frequency conductivity 
image spectra with high spatial resolution and sensitivity using MR scanner within the 
range of DC to the Larmor frequency since significantly distinguishable signal changes in 
biological tissues depending on the applied frequency range. 

We have a plan to support the clinically applicable combined conductivity imaging 
method based on MR techniques. We may apply this method to detect the destroyed cell 
membrane or the modification of cell or tissue structure due to necrosis or apoptosis ini- 
tiated by inflammatory response inside the body. Since RF ablation and cryoablation for 
cancer treatments are accompanied by the destruction of cell membranes, it will be a good 
method to diagnose ablated lesions based on physiological status of tissue to improve the 
safety and predict the local recurrence after ablation. 
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Conclusions 

Cross-sectional conductivity imaging methods for high spatial resolution and sensitivity 
inside the human body has been actively investigated. MREIT and MREPT techniques 
show different electrical properties at low frequency (below 1 kHz) as MREIT uses an 
externally injecting current while MREPT uses the Larmor frequency of an MRI scanner. 
Recently, the dual-frequency conductivity imaging from a single MR scan simultaneously 
was proposed based on combination of MREPT and MREIT. Even though it produces 
conductivity spectral information, MREPT and MREIT have commonly suffered from 
weak signals and noise amplification since the procedures to reconstruct the conduc- 
tivity in MREIT and MREPT need to differentiate the measured phase signals. We 
suggested the optimization method to find weighting factors according to echo signals for 
MREPT and MREIT using the ICNE multi-echo pulse sequence, which minimized the 
noise artifacts in the measured phase data. The noise standard deviations were precisely 
analyzed for the optimally weighted magnetic flux density and the phase term of positive- 
rotating magnetic field. We applied the denoising method based on the reaction-diffusion 
equation with the estimated noise standard deviations to enhance the quality of dual- 
frequency conductivity images. A real phantom experiment was performed to validate 
the proposed method using common measured data to reconstruct the dual-frequency 
conductivity distributions for MREPT and MREIT. 
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